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Abstract. We review the random matrix description of electron transport 
through open quantum dots, subject to time-dependent perturbations. All 
characteristics of the current linear in the bias can be expressed in terms of 
the scattering matrix, calculated for a time-dependent Hamiltonian. Assuming 
that the Hamiltonian belongs to a Gaussian ensemble of random matrices, we 
investigate various statistical properties of the direct current in the ensemble. 
Particularly, even at zero bias the time-dependent perturbation induces current, 
called photovoltaic current. We discuss dependence of the photovoltaic current 
and its noise on the frequency and the strength of the perturbation. We also 
describe the effect of time-dependent perturbation on the weak localization 
correction to the conductance and on conductance fluctuations. 
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1. Introduction 

Quantum dot is a small disordered or irregularly shaped conductor, connected to 
leads P, see Fig. 1. Exact values of the conductance of a quantum dot are determined 
by electron wave functions in the system and are hard to calculate exactly for arbitrary 
configurations of the dot. Moreover, the conductance changes significantly even for 
tiny changes in the position of impurities or the boundary of the dot. Due to extreme 
sensitivity of the conductance on many parameters, the statistical description of the 
conductance is more appropriate [SEIEEIEIIZIIHIEI • The random fluctuations 
of the conductance from sample to sample of non-interacting systems are universal. 
The universality jllL [T^ means that the conductance statistics can be described by 
universal functions, which are independent from the shape of the dot or the details 
of the disordered potential. Particularly, the variance of the conductance varg is of 
the order Gq and is nearly independent from the sample geometry (Go = e'^/irh is the 
quantum of conductance for spin degenerate electrons). The other universal quantity 
is the weak localization correction to the conductance, defined as the difference of 
average values of the conductance over orthogonal (zero magnetic field) and unitary 
(strong magnetic field) ensembles. The weak localization correction to the conductance 
is also of the order ofGoHSHHIITSllTniinillllini. 

A common description of electron transport through quantum dots is based on the 
Landauer formalism (20) I21L 1221 1^ , when the transport characteristics of the system 
are described in terms of the scattering amplitudes between different conducting 
channels in the leads. There are several approaches for statistical description of 
electron transport. One approach is based on a diagram technique developed for 
disordered bulk metals |23i when the scattering amplitudes are represented in terms 
of electron Green functions H ESI EEl 123 ■ 

Alternative approaches are based on the description of the system by random 
matrices, when either an exact scattering matrix is replaced by a random unitary 
matrix, or an exact Hamiltonian is replaced by a random Hcrmitian matrix. In the 
first case, the unitary matrix is taken from Dyson's circular ensemble of uniformly 
distributed random matrices |28[ll8[[Tn] . In the Hamiltonian approach, the Hermitian 
matrix belongs to an ensemble of random matrices '29' with the Gaussian distribution 
of its matrix elements [301 12] ■ The equivalence for statistical description of electron 
transport by both random matrix approaches was shown in |31l 1321 133| . 

Although the random matrix approach is not based on microscopic description 
of electron system, their correspondence to microscopic problem has been proven for 
disordered metal grains |34[ I35L . The validity of such random matrix description 
of chaotic ballistic systems was addressed in |^ ■ 

We imply the following realization of the system, see Fig. 1 a). Negative voltages 
applied to the gates (black areas) confine electrons to a small region (light gray), 
forming a quantum dot. Electrons in the dot are connected to the electron reservoirs 
by narrow leads. Electric current that flows through the dot can be measured as 
a function of the voltage bias V between the reservoirs and the amplitudes of ac 
gate voltages Vi^2{t)- Particularly, the current linear in bias V is determined by the 
conductance of the dot. Changing magnetic field or shape of the dot one can obtain 
different realizations of the quantum dot and experimentally study statistics of the 
quantum corrections to the conductance. 

The quantum corrections to the conductance are commonly characterized by the 
weak localization and the variance of conductance fluctuations. As any other quantum 



Quantum chaotic scattering in time- dependent external fields 



3 



interference phenomena, they are very sensitive to inelastic processes, commonly 
referred to as dephasing [3^1 • A phenomenological description of the effect of dephasing 
on electron transport through open quantum dots was developed in Refs. 001 ■ 
The dephasing rate due to electron-electron interaction in quantum dots was estimated 
in 1331 121 • Another possible source of dephasing is a time-dependent perturbation, 
such as a microwave radiation or periodic deformation. In this case the Hamiltonian 
of the system can be considered as a time-dependent random matrix _45 , 46 , 47 , and 
all transport quantities can be calculated as a function of various parameters (e.g. 
strength and frequency) of the time-dependent perturbation. The scattering matrix 
description of the system subject to time-dependent perturbation was developed in 
energy representation by Biittiker, Thomas and Pretrein in Refs. j^HlEH. In this case 
the scattering matrix describes processes when electron scattering between different 
channels in the leads is accompanied by the change of electron energy. Alternatively, 
the analysis of the effect of time-dependent perturbation on the conductance can be 
carried out in time representation, see Refs. 03 0^1 ^1 1^01 In general, both 
the weak localization correction to the conductance and the variance of conductance 
fluctuations are suppressed by time-dependent perturbation. The suppression of 
the quantum corrections to the conductance by microwave radiation was observed 
experimentally in 

Time-dependent perturbation of quantum dots not only suppresses quantum 
corrections to the conductance, but also produces electric current through the system 
even at zero bias. This effect is related to the charge pumping, which occurs in systems 
with large tunnel barriers |52l 1541 EHl ■ If the conductance of the system is very 
small, the electric current is quantized in units of ea;/27r, where 27r/w is the period of 
the pump. At finite conductance, a countercurrent reduces the pumped current and 
thus violate the quantization of electric current |56| . For an open quantum dot, the 
countercurrent nearly compensates the pumped current and the current is no longer 
quantized. 

In the low frequency limit, the magnitude of the pumped current is determined 
entirely by the evolution of the system in the parameter space, see Fig. 1 b), under 
time dependent perturbations |^ EHl EHl EDI • As frequency increases, the parametric 
description becomes insufficient and requires full analysis of electron dynamics in time- 
dependent fields [fill I47| . The analysis of how the adiabatic description breaks down 
at finite frequency can be also found in Refs. [^21 We note that the charge 

pumping through an open quantum dot is a manifestation of the photovoltaic effect, 
which occurs in systems without inversion center |64| . The photovoltaic effect was 
previously considered by Falko and Khmelnitskii |HS1 in mesoscopic microjunctions 
and by Kravtsov, Aronov and Yudson |^ in normal metal rings. 

It turns out [^3 that the photovoltaic current is sensitive to the actual electron 
distribution function in the dot. Time-dependent perturbations may broaden the 
distribution function, resulting in heating. This broadening of the electron distribution 
occurs as a result the electron diffusion in the energy space. The effect of 
time-dependent perturbations on electron distribution function becomes even more 
interesting in closed systems, when the energy diffusion acquires quantum interference 
corrections. The latter leads to a dynamic localization of the electrons in energy 
space El|7ni. 

Photovoltaic current fluctuates not only with respect to different realizations 
of the quantum dot, but also for a given realization due to quantum and thermal 
fluctuations. Such fluctuations are called current noise and are described by the 
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Figure 1. a) Schematic picture of the experimental setup. Light gray color shows 
the region available for free electron motion, while the dark gray color shows the 
region forbidden for electrons due to electrostatic repulsion from the gates (shown 
in black) with applied negative voltages. A finite bias V is applied between the 
upper and lower (in the text refereed to as left and right) reservoirs. Oscillating 
voltages Vi^2{t) applied to the gates produce time-dependent perturbation of the 
electron system, b) Contour plot represents time evolution of gate voltages Vi_2(t)- 



fluctuations of the charge transported through the dot in a certain number of 
perturbation cycles. The statistics of such charge fluctuations was studied in Refs. 
| 71l I72L 1 731 1741 175L 17^ for temperatures T and pumping frequencies u; much smaller 
than the inverse dwell time 7csc (escape rate) of electrons from the quantum dot. 
Particularly, Refs. O E21 addressed the full counting statistics at temperatures 
T <C w (we use h = 1 and the Boltsmann constant fee = !)• The mean square 
charge fluctuations for uj,T <^ 7esc (but for arbitrary relation between oj and T) were 
considered in Ref . 76 . The variance of the photovoltaic current for arbitrary relation 
between the temperature T, the frequency w, the escape rate 7csc and the strength of 
the perturbation was calculated in Ref. J77,. 

Experiments [751 [ZHl IHOI were performed to detect the photovoltaic current in 
various mesoscopic systems [7S1[7H], including open quantum dots W in the adiabatic 
regime. The observed magnetic field symmetry and the amplitude of the current 
indicate that the measured current was likely related to the ac rectification PT1I5^I5S| . 
A more detailed analysis of the zero-bias current in different regimes of microwave 
radiation shows that in some instances the photovoltaic current, and not the 
rectification current, was observed |H21E1- 

In this review we focus on the random matrix description of electron transport 
through open quantum dots in the limit of the large number of open channels A^ch 
connecting the dot to the leads. This condition allows us to neglect the electron- 
electron interaction that gives corrections of the l/A^jj order, see Ref. The 
same condition permits the use of a diagrammatic technique, similar to that described 
in [23, to calculate ensemble averaging. We assume that the electron dynamics in 
the dot is fully chaotic and disregard classical fluctuations of the conductance j26j . 
We emphasize that the random matrix description is applicable for sufficiently small 
quantum dots, when the Thouless energy = l/rcross is much greater than all 
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other energy scales of the problem, such as the frequency lu of the perturbation or 
the temperature T (tci-oss is the electron crossing time of the dot). Larger systems 
{E^ <C w, T) can be treated by methods developed for bulk conductors '3^, see 
e.g. [SU]. We note that the derivation of the results will be performed within the 
Hamiltonian formalism, following Refs. 021 1^ i but the same results were derived 
within scattering matrix formalism in |47j . 



2. Scattering Matrix Formulation of Transport through Open Quantum 
Dots 

2.1. Model 

The Hamiltonian of the system is 

Ti, — Tid + Ti-id + Til- (1) 
We choose the basis for electron wave functions in the dot, so that the coupling of 
states in the dot to states in the leads can be written as 



T, if n = a < iVch, 



^ — 0, otherwise. 

Here ipn and ipa{k) are the annihilation operators of electrons in the dot and the 
leads, respectively. Index n enumerates electron states in the dot: n = 1, . . . M, with 
M —i oo. Index a labels channels in the leads, with 1 < a < A'j for the Ni channels 
in the left lead and with A^i + 1 < a < iVch for the N,- channels in the right lead, 
■^ch = N\ + N,-. Coupling constants T are defined below in Eq. H43|l . The Hamiltonian 
for electron states in the leads near the Fermi surface can be linearized: 

Fi = «F^fcV'l(fc)V'a(fc), (3) 

where the continuous variable k denotes electron momenta in the leads, vf = (2tti')~^ 
is the Fermi velocity, and v is the density of states per channel at the Fermi surface. 

Finally, Tid is the Hamiltonian of the electrons in the dot, determined by the 
M X M matrix H and the electrostatic energy of N electrons: 



rp + E,N^. (4) 



Matrix H describes the time independent part of the electron Hamiltonian, and the 
time-dependent component of the Hamiltonian is represented in terms of the traceless 
matrices Vi and the diagonal matrix IZi. In the setup shown in Fig. 1, the time- 
dependent perturbation is generated by the gate voltages Vi^2{t). The perturbation is 
linear in small amplitude of oscillating voltages Vi,2(i), and the time evolution of the 
perturbation characterized by the dimensionless functions (pi{t) oc Vi{t). The second 
term in Eq. represents the largest in 1 /M contribution from the electron-electron 
interaction with Ec being the charging energy of the dot, and N = being 
the operator of the electron number in the dot. The status of this approximation 
was discussed in detail in Ref. j85j . For an open quantum dot with the large number 
of open channels N^h ^ 1 the interaction term can be treated within mean field 
approximation, and the Hamiltonian Ti.d in Eq. can be further simplified: 



tA, eVd{t)=J2z^^^it) + ^E,{N). (5) 
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Here we introduced the electric potential Vd(t) linear in the quantum mechanical 
average (N) of the electron number N in the dot. Corrections to this mean field 
treatment were calculated in Refs. [MII5^I57| . 

To determine the electric potential Vd(t) in the dot, we have to define the quantum 
mechanical average {N) of the electron number iV in the dot. In each particular 
moment of time the electron number (N) is not constant and its time evolution is 
described by the discontinuity equation eN{t) = Ii{t) + I\{t). We estimate Vd{t) to 
the lowest order in 1/iVch ^ 1, and use NiGq (Ni-Gq) for the conductance of the left 
(right) contact, Gq — c^/tt is the quantum conductance. Then, {N{t)) satisfies the 
following equation 

^ = -7e.(iVW) + ^ (1^ - V,it)) + ^ - VM) . (6) 
at TT n 

The first term in Eq. ||HJ) is a diffusion term, describing the electron escape from the 

dot with rate 7csc, where is the mean time for an electron to escape the dot 

through one of the leads; below we define 7esc in terms of microscopic parameters of 

the system. The last two terms in Eq. © represent electron flux from the dot due 

to the voltage difference V|.(r) — Vd{t) across the contact of the left (right) reservoir 

and the dot. A discussion of the charge dynamics in quantum dots can be also found 

in [HHl. 

Combining Eq. © with the expression for Vd{t) from Eq. ©, we obtain 



The characteristic energy scale governing the dynamics of the charge is EcNch/^n cx 
GoNch/Cd, Gd is the dot capacitance. Usually, this scale is of the order of the Thouless 
energy Et and significantly exceeds electron escape rate 7osc- Therefore, we consider 
the limit, when both 7csc and the frequency of the external field to are much smaller 
than Ec, and use the following equation for the electrostatic potential of the dot 

Vd{t) ^Vd= ' '^^^ ^ ^ (8) 

We conclude that the time-dependent perturbation Eq. |0J can be chosen traceless, 
Zi = 0, and the electric potential in the middle dot is determined by the potentials of 
the left and right reservoirs. 

2.2. Electric Current 

The current through the dot is given in terms of the scattering matrices S{t, t') by the 
following expression 

jdtidt2tv[k[s{tM)I{ti^t2)S\t2.t)- f{+Q)\Y (9) 

The derivation of Eq. can be found in j49l IHll Wf\ , see also [Appendix A| Here 
(/) stands for the quantum mechanical and thermodynamic averages of the current 
operator (no ensemble averaging!) and f{t) represents the electron distribution 
function in the leads in time representation. We consider the case when electrons 
in the leads are in thermal equilibrium at temperature T, but the different voltages V\ 
and Vr are applied to the left and right electron reservoirs. Then, the matrix f{t) is 
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diagonal faa{t) = fi{r){t), if channel a belongs to the left (right) lead. The function 
f\(r){t) is the Fourier transform of the Fermi-Dirac distribution function: 

+00 

/,,„(.) ^ .""'-/(r); fir) ^ I 1^ e'- - 1} ^ (10) 

— 00 

Here the traceless diagonal matrix A is introduced 

if 1 < a < M; 

= { (11) 

^, ifiVi<a<iVch, 

JVch 

and the scattering matrix S{t, t') 

SMt,t') = e-^<^(*-*') [S^pS{t - t') - 27:z,yWLG^„,{t,t')Wmp] , (12) 
is defined in terms of the Green function G^„^{t, t') that satisfies the following equation 

\i^-H-Y,y^V^ii) + ^^^^^^] &''Kt.t') = 5{t-t'), (13) 



where the matrices H , Vi and W were introduced earlier, see Eqs. |5J) and lO- The 
diagonal component eVd of the electron Hamiltonian in the dot is removed from the 
expression for the electron Green function G^„^{t,t') by the gauge transformation, 
represented by the exponential factor in Eq. (|12|l . 

To the linear order in voltage across the dot V = Vi — Vr, the dc electric current 
has the form 

(/)=/ph + gF. (14) 

The first term represents the photovoltaic current, which flows through the dot even 
at zero bias. The second term is linear in voltage V with factor g being the dc 
conductance of the dot in the presence of time-dependent perturbations Vi. The 
linear in V contribution to the current in general may come from two sources: i) the 
non-equilibrium distribution of electrons in the leads and ii) change in the photovoltaic 
current /ph due to change in the configuration of the electron wave functions when 
the bias is apphed. Due to the electro- neutrality condition Eq. {71), the voltages Vi, 
Vd and Vd enter only as exponential factors to the expression for the electric current 
Eq. ©, and do not actually affect the structure of electron wave functions in the dot. 
Therefore, only the non-equilibrium current contributes to the linear in V term. 
The dc conductance g of the dot can be represented in the form 

To +00 

g = gci + Go fdtidt2F{ti~t2)tT[s{t,ti)AS\t2,t)AY 5ci = Go^^. (15) 







Here g^ is the classical conductance of the dot, Tq is the observation time, Gq ~ e /ttH 
is the quantum conductance for doubly degenerate electrons in spin states and F{x) 
is the Fourier transform of the derivative of electron distribution function: 

Fit) = (16) 
^ ' sinhTrTi ^ ' 
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The expression for the photovohaic current /ph can be obtained from Eq. @ by 
taking Vi = V^. Using the Wigner transform for the scattering matrix 

Sit,t') = ls^ie)e^<'-''^^, (17) 

we write 

To 



For slow perturbations ipi with frequencies uji smaller than temperature T or the 
inverse eigenvalues of the time delay matrix | (5^ 

n,{e,t) = [dA{£)]Sl{e), (19) 
we can expand the scattering matrices in Eq. IjlSj) in r and obtain 

7dt f de 1 i 7 fdSt(e) . , , ^ , , dS} (e) \ ] 



To 



Here, the scattering matrix in the Wigner representation is a function of the 
perturbation itself and its time derivatives: St[£) — S{e, (pi{t), ^Pi{t), . . .) because the 
Green function G^^^ (e, t) is a solution of the equation: 



rG(e,t) - i -iTTiyWW^;G^^'>{e,t)^ 

1 d'^'P.it) (yd 
2{2iYk\ df" [ " de>' ' ' de'' 



fe = i 

with {A; B} = AB + BA. In the adiabatic approximation the derivatives d'^ ipi{t) / dt'^ 
can be neglected and the scattering matrix is determined by parameters (pi{t), 

, . /^ tr{AIm7^.(£,(^)} ^^ dSie^ 
-^ph = TTT 9 d(Pi / — ,9 ^»(£,<^) = — '5T(£,(p). (22) 



/ph T^j-^^j cosh2£/2T 2^' 

The integral in Eq. H22|l runs over the loop in the parameter space (pi, and Tp is the time 
for a system to complete this loop. Particularly, for the perturbation characterized by 
two parameters 

tpi{t) = Xi cos(wt), (p2{t) ~ X2 cos(ijji + 0) (23) 
the photovoltaic current is given by |57| 



eoj 



dS dS^ 



where the integral runs over the inner part of the ellipse, defined by Eq. (|23|l . At finite 
frequencies, but still uji <C T, Eq. H20|l is still applicable and can be rewritten in the 
form similar to Eq. H22|l . if the parameter space (pi is extended to the phase space, 
containing time derivatives of fi{t) as well. 

With the help of the equations of motion Eq. l(T^ . the expression for the 
photovoltaic current can be rewritten in terms of the Green functions G^'^{t,t'): 

To 

/ph = 2ez^z/ J Jdtidt2Fl'^{ti,t2)Y,^r[w^G'^{t,ti)V,G^{t2,t)WAY (25) 
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which is more convenient in some calculations. Here the function 

i^P^tl - t2) = [(p.(ti) - ip,it2)] fiti - t2), (26) 

takes into account the probability of electron transitions due to the perturbation 
V(pi{t) for the equilibrium electron distribution in the dot f{t). We notice, that 
taking higher order terms in Vi in G^'^ results in a new electron distribution function 
in the dot: 

V,Fl'^{t~t') G^(t,ii)l/,i^f (ti -t2)G^(i2,i')- (27) 
The effective electron distribution function has a shape different from the Fermi 
distribution function, see Sec. [S] 

We note that due to trA = the expression for the conductance cannot 
be represented in terms of modified distribution function. As a result, see j46| . 
conductance fluctuations are characterized by the electron temperature in the 
reservoirs rather than by the electron temperature in the dot. This statement was 
further investigated in j51| , where the effect of time-dependent perturbation on two 
possible definitions of the conductance was studied. It was shown, that the Landauer 
conductance, defined as the linear response to the bias between the reservoirs and 
given by Eq. H15|l , is indeed characterized by the electron distribution function in the 
leads. In other geometries one can measure the linear response of electric current 
to the internal perturbation of the mesoscopic system by the dc electric field. Such 
response, called the Kubo conductance, is sensitive to the actual distribution hmction 
of electrons in the mesoscopic system. 



2.3. Current Noise 

The current correlation function S represents fluctuations of the charge Q — /J"° I{t)dt 
transported through the dot over the observation time interval Tq 

s = M^M! = r amm) - {m){m)) (28) 

Expression for S in terms of the scattering matrices S can be derived in a similar way 
to the derivation of Eq. @ for current, and is outlined in [Appendix B| For arbitrary 
distribution function fa(3{t) = Sa/3fa(t) in the leads the current correlation function 
S has the form {5t,t' = 6{t - t')): 

8=1^° dtdt' J dtidt2dt[dt'2 tr | (sHt2,t)AS{t, t[) - kSt,,tStM^'^ fit[ - t'^) 

X (sHt'^,t')kS{t\t^)~k5e^^t,6t>M)i^5uM- hii-t2))]. (29) 

Below we consider the case of zero bias across the dot, so that fa[t) = f{t), see 
Eq. H1U|) . We also assume that the temperature of the system T is finite and Ttq ^ 1. 
(The limit T = has some interesting properties and was discussed in Refs. [T^ITT] '). 
Then, S can be divided into two parts: 

S = Snj + Sp. (30) 

Here, the second term Sp is chosen in such a way, that in the absence of time dependent 
perturbations this term vanishes [^p = 0, see Eq. (|32|l below], and only the first 
term remains. The first term describes the current noise due to thermal fluctuations 
of electrons in the leads at temperature T and is known as the Nyquist- Johnson 
noise [5011^ . 
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The Nyquist-Johnson component of the noise can be written as 

To 

Sn, = 2g,{r^ j ^ j dtMit, ^ t')f{t' - t2)tT [AS{t,t,)ASHt,,t)} , (31) 



where f{t) — 6{t) — f{t). The first term in Eq. (|31|l represents the noise of a classical 
resistor with resistance l/gci- The second term in Eq. H31|l describes the contribution 
to the current noise from the quantum mechanical corrections to the conductivity, cf. 
Eq. IjlSfl . In the absence of time-dependent perturbations, the second term represents 
the quantum correction to the conductance, so that the noise correlator has the form 
Snj = ZgT, where g is the sample-specific conductance of the dot, see Eq. (|15ll . 

The external field changes the conductance of the dot, see Sec.0] Consequently, 
we can expect that the Nyquist-Johnson contribution to the current noise is also 
modified due to the external field. In particular, the ensemble average 5nj and 
fluctuations of S'nj with respect to different dot realizations are suppressed by time- 
dependent perturbation. 

The second term, 5'p, in Eq. H30(l represents the noise of the photovoltaic current 
Eq. I)18|l and has the following form in terms of the scattering matrix S{t,t'): 

Sp = e2 r ^ / dtidt2dt[dt'Jiti - t2)f{t[ - t'a) 



X tv{s{t'^,t)kS^{t,t,)Sit2,t')AS^it',t[) - AHt'^^tSt,tA,,t'St',t'^}. (32) 

To discuss the noise of the photovoltaic current in more detail, we consider the 
adiabatic limit, when the eigenvalues of the time-delay matrix Eq. 119(1 are shorter than 
both 1/T and l/tOi, {uji is the frequency of external perturbation Vi). The ensemble 
average value of Sp for arbitrary strength and frequency of the perturbations was 
investigated in [771 IHl] and is briefly discussed in the end of Sec. 

In the adiabatic limit only electrons close to the Fermi energy contribute to the 
current. Thus, we can neglect energy dependence of the scattering matrix S{e, t) in the 
Wigner representation Eq. ((17(1 and substitute S{t,t') = St{e = 0)St,t' into Eq. 1(32(1 : 



dtdt 



Sp ^e'j^ = tr {A^ ~ sU0)ASMStmSA0)} fit - t')f{t - t'). (33) 

We observe that the temporal correlations in the current survive on time scales 
comparable with the observation time Tq at low temperatures Ttq ^1. In this case 
the full counting statistics is non-trivial and higher moments of the current should be 
investigated, see Refs. for more detail. At finite temperature T, the temporal 

correlations of the current are suppressed on time scale of the order of 1 /T, see Eq. ((l()(l . 
and the counting statistics of the current becomes Gaussian and is described by the 
average value of the current /ph, Eq. 1(18(1 and its noise 5p, Eq. 1(32(1 . 

Within a bilinear response to the perturbation Eq. ((23(1 we obtain the following 
expression for the noise 

Sp = e^Fn{T,u) [JCiiXl + JC22XI 2COS0 . (34) 

Here coefficients /Cy are given by 

/C,, =tr{[A;7?:,][7?:j;A]} (35) 
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and the function Fa{T,u}) represents the probabihty of the absorption or emission of 
a perturbation quantum with energy uj: 

^ [/(£ + ^u;)f{e - i^) - 2/(e)/» + f{e - ^hi.)f{e + ^Tico) 

At low temperatures T <C i^, but stiU T ^ I/tq, Fn = lo/2tt. As T increases, Fn 
decreases Fn = w^/T. 

Above we discussed the current noise in the situation when the bias across 
the dot is zero. When a finite bias is apphed, the noise acquires dependent on 
the bias contribution called shot noise. It was shown |M] that the shot noise 
originate only due to quantum corrections to electron transport, while the classical 
contribution to the transport does not lead to the shot noise. Therefore, one can 
expect that a time-dependent perturbation suppresses shot noise along with any other 
quantum interference characteristics of electron transport. Another interesting effect 
of microwave radiation on the shot noise of open quantum dots was found by Lamacraft 
in Ref. [M]. This effect results in cusps of the noise power when the bias eV is a 
multiple of microwave frequency oj: eV = nu with integer n. 



3. Ensemble of Open Quantum Dots 

The exact form of the Hamiltonian Eq. Q for quantum dots depends on many 
microscopic parameters of the system, such as the shape of the dot, position of 
impurities and is usually too complicated for analysis. However, for many purposes the 
interesting question is what the statistical properties of transport coefficients through 
a quantum dot, rather than the corresponding values for each particular sample. To 
describe statistical properties of quantum dots, a random matrix theory turns out to 
be a productive tool. The random matrix description of quantum dots is based on 
the assumption, that the Hamiltonian of the dot, Eq. is determined by M x M 
matrices H and Vi with H being a random realization of a hermitian matrix from the 
Gaussian ensemble The matrix elements Hnm{^) of matrices from this ensemble 
in the presence of magnetic flux <I> through the dot are described by the following 
correlators 

Hnm(.^l)Hl,^,{^2) = ^ [F($l - <^2)5nn'5mrn' + + $2 )<5™„' (5„™' ] , (37) 

Here (. . .) stands for the ensemble averaging, and 5i is the mean level spacing of 
eigenvalues of H . For small A<i> function F(A<i>) can be estimated as F(A$) — 
1 — K(A$/<i>q) , where k is a non- universal, sample-specific constant of the order 
of unity, and $q = c/e is the flux quantum [HHIIZI- At $1^2 — 0, the matrix i?(0) 
belongs to the orthogonal ensemble. As A$ increases, L(A$) vanishes, and H{^) 
becomes a matrix from the unitary ensemble when the second tern in Eq. 137|l is equal 
to zero. The microscopic justification of the random matrix description Eq. 1)37(1 can 
be found in I35L 15^ for disordered systems and in EHl in ballistic chaotic 
systems. 

Matrices Vi can also considered as Hermitian random matrices. Below we 
disregard the fluctuations of the matrices Vi, and assume that Vi are real symmetric 
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matrices, belonging to an orthogonal ensemble. In this case, we characterize 
perturbations Vi by parameters 

C., = ^trl/,l/,. (38) 

We remind that trV^ = 0, sec Eq. (|7|). The parameters have the meaning of the 
level velocities which characterizing the evolution of an energy level Eni^) under the 
external perturbation J2i Vi(pi{t), [HIHITT] : 



TT difi dipj dipi dipj 

Parameters are also related to the transition rates of electrons under perturbation 
Vi. Indeed, the transition rate 711 due to perturbation Vi is determined by the Fermi 
golden rule: 



= V 2n\V,..n^\H{e^ - ± c.) ^ %^ ^ ^. (40) 

Ol TT 

m 

The first equality sign follows from the Fermi golden rule, the second sign represents 
an estimate of the characteristic value of the matrix elements |Vi;„mP and the density 
of states l/i5i, the last equation is the definition of Cn, cf. Eq. 1)38(1 . If perturbations 
induce uniform electric fields Ei in a quantum dot with typical length L, parameters 
Cij can be estimated as Cy ~ e'^EiEjL'^ / ETh, where i?Th ~ MSi is the Thouless 
energy. 

Below we show that all statistical transport characteristics of quantum dots 
in the presence of time-dependent perturbations are functions of parameters Cij. 
Thus, even though Cij are free parameters, measurements of several transport 
characteristics ^1 EHl allow one to eliminate the uncertainty of Cij . 

For calculations of different correlation functions of transport parameters over 
the ensemble of random Hamiltonians -ff , we use a diagrammatic technique, similar 
to one developed for disordered metals, see In this Section, we briefiy discuss the 
basic elements of this diagrammatic technique. 

First, we calculate the ensemble averaged Green function G^'^{e) in the absence 
of time-dependent perturbations. The diagram equation in Fig. |21 reduces to the 
following algebraic equation for the electron self-energy I](e) ~ {M61/tt'^) tr G^{e): 

^ _ MSl 1 M6f 1 1 

7r2 e_S(e)+iO s - ^{e) + lO e ^ ^(e) + iMSi/ir' ^ ' 



Solving Eq. (|41|l . we find the ensemble average Green function G^{s) = {G^{e))* for 
£ <ti MSl in the form 

Nch + i2ne/Si 
+ 4M ' ^ch < n < M; 



GS„M ^ -^S^njj- \ ^ (42) 
In derivation of Eqs. ((41|l and ((42(1 . we used the following values for factors T in Eq. ^ 



V TT^V 

This choice of T corresponds to a dot connected to the leads by reflectionless contacts, 
when the ensemble averaged scattering matrix Sap is zero and S belongs to circular 
ensemble, see [7]. 
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+ iTTuWW^ " 



Figure 2. Diagrams for the ensemble averaged electron Green function in 
the dot. The first line of the figure introduces the bare Green function [e + 
i-KvWW^]~^ and correlation function of the matrix elements of the Hamiltonian 
H. The second line represents the Dyson-type equation for the ensemble averaged 
Green function G^„(e). The third line introduces the first two terms of the self- 
energy S(e), which is diagonal in index of electron states in the dot. The second 
term as well as all other terms which contain intersections of dashed lines are 
small in parameter 1/M. 



We also introduce two other elements of the diagram technique used in 
calculations of statistical properties of electron transport in the presence of time- 
dependent perturbations. One element is called the diffuson t2, t) and is defined 

by 

4M2A2 

^-D^.^^A'-^^'-^^A-tt). (44) 



2 ' 2 

see Fig. 121 a). The other element is called the Cooperon C(Ti,T2,t) and is defined as 

[Gf.rn-,^Srt,tt)Gi^,^^{t-,,t^)]^^^ = ^^5{t+ + +t^ t\ - tr) 

xC$,+*, (^t+-tr,t2+-<2",^^^) , (45) 

see Fig. 121 b). These two elements represent the ensemble average product of the 
advanced and retarded electron Green functions in the dot, divided by the product 
G^n;^r Gmm;*i ^^^.^^ G^^.^^ , so-callcd the amputated average. The diffuson and 
the Cooperon are given by the following expressions: 

2?A*(ii,t2,T) =0(ti-i2)exp JrA<i>{T,t)dt^ ; (46) 

CA*(ri,T2,t) =e(Ti-r2)exp ^-i JrA*(T,i)dTj . (47) 
Here we use the notation 

rA*(r, t) = 7e.c + 7(A$) + ^^i^^ t)G^3^J{T, t), (48) 

7csc = ^^^; 7(A$) = -[l-L]; Mr,t)=^;{t + T/2)-Mt-r/2).m 

Ztt tt 
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Figure 3. The Dyson- type equations for the difTuson, Eq. 1441 . and the 
Cooperon, Eq. 1451 . are shown in the first and the second fines respectively. The 
subscript "amp" in Eqs. 1441 and I45i emphasizes that the four Green functions at 
the terminals of the diffuson and the Cooperon are omitted. All Green functions 
in these diagrams and diagrams in figures below are ensemble average Green 
functions, introduced in the second line of Fig. 1^ 



The first term in Eq. (|48|l is the electron escape rate from the dot and 7a* is the 
electron dephasing rate due to the difference in magnetic flux A^. The last term 
in Eq. 148|l describes the effect of time-dependent field on the correlation functions 
Eqs. (jHll and of electron propagators. Equations and (|T7|) are the 

building blocks of the diagrams, which are studied below for different correlation 
functions of transport characteristics of open quantum dots. 



4. Effect of time-dependent perturbations on the conductance 
4-1- Weak localization correction 

Weak localization correction to the conductance of a quantum dot is given by the 
ensemble average of the second term in Eq. H15|l . For the unitary ensemble weak 
localization correction (23 is small as f/ci/^ch 1 beyond the accuracy 

of our calculations. In the orthogonal ensemble the weak localization correction is 
gc\/Nch ~ 1 [Z|- We define the weak localization correction to the conductance as the 
difference between the averaged values of the conductance over orthogonal ($ = 0) 
and unitary (<i> ^ $q) ensembles: 



Afiwi = (g)<i>=o - (g)<i>»<i.q- (50) 

In this subsection we describe the effect of time-dependent field on the weak 
localization correction Eq. H5U|) . 

The weak localization correction is given by the diagram in Fig. ^ and can be 
calculated from the following expression j45| : 

27r/0 oo 

A5wi = Ag^^ J ^ Jj^^^drCir, -r, t), A^i? = -Go^- (51) 



This equation gives the universal description of the effect of the time-dependent 
fields on the weak localization correction. Below we will discuss different asymptotic 
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Figure 4. The diagram for the calculation of the weak localization correction 
to the conductance. 



regimes for the case when the perturbation is described by only one harmonic function 
ipi{t) = coaut and Cn — C\. 

In the absence of the time dependent perturbation Ci = 0, one obtains Ag^\ = 
^g'^i mUHI' For weak external field Ci ^ 7csc we find 

where 7csc is defined in Eq. H49|l . In this regime the correction is quadratic in the 
frequency of slowly oscillating field, similarly to the result for bulk metal system at uj 
smaller than the dephasing rate l/r^. However, the frequency dependence saturates 
at large frequency. It is different from the result for bulk systems where a 

characteristic spatial scale shrinks as y/D/oj with D being the diffusion coefficient, 
whereas in a quantum dot this scale is determined by the size of the dot L. The 
random matrix description breaks down at ^ D/L^ = -Etii, £'Th is the Thouless 
energy. 

In the opposite limit of strong external field Ci 3> 7csc we consider separately the 
cases of fast, uj ^ ^csc, and slow, uj <C 7gsc oscillations. In the first case we have 



Agwi /7csc .^„s 

(53) 



The 1 / \/C\ power dependence of the quantum correction is similar to that for the bulk 
system. Contrary to the bulk systems, the result does not depend on the frequency oj 
for the reason mentioned above. In the case of slow field uj <C 7esc, but still Ciw^ » 7^^^ 
(strong field) the weak localization correction to the conductance is 



Affwi r(i/6) 3 

The power and frequency dependence of Ag^i is again different from that in bulk 
disordered metals, r(a;) is the F— function. 

4-2. Conductance fluctuations 

Next we consider the fiuctuations of the conductance g over ensemble of random 
Hamiltonians H. We notice that the fiuctuations in g originate only from the second 
term in Eq. H15|l . since in the model of fully chaotic quantum dots with open channels 
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the classical conductance does not fluctuate. We have for 5g = g — gd the following 
correlation function, which can be derived j46| from the diagrams shown in Fig. [S] 

27t/lu oo oc 

^9^Jm:=^ll. I ^^^^ fdrF^r) f dO [K+ {t,t' ,T,e) + R- {t,t' ,t,0)] , (55) 

ch J J 

r/2 

where the functions K^{t, t' , t, 9) are given by 

^-(M',r,.)^p(i±^,^±4±^-M'-i)p(^,^4^-M-t'), (56) 

K~{t,t\T,e)^c(^t-t' + e-^,t-t'-9 + ^,^-±^^ + ^ 

xc(t'-t + e+^-,t'-t-e-^-,'-±^^^y (57) 

The two terms in Eq. H55(l have different properties with respect to the magnetic 
flux $ through the dot. Although at $ = both terms survive, at finite magnetic 
field, |$i,2| ^ ^q, only one of them remains: for $2 ~ $1 the second term vanishes, 
and for $2 ~ —^1 the first term vanishes. The values of the conductance correlation 
function at $2 — ±$1 characterize the symmetry of the conductance with respect 
to magnetic field inversion. If conductance is symmetric, Sg<i6g^ — Sg^dg-<i. This 
equation is indeed valid in the absence of time-dependent perturbations [Hnil221- As 
was shown in Refs. |46II50) . time-dependent perturbations may suppress the symmetry 
of the conductance with respect to magnetic field inversion, see e.g. Eq. (|61|) below. 

In the presence of a single harmonic perturbation at frequency uj and with strength 
Cii = Ci, we can write the conductance correlation function in the form 

■%<3.f£L,i.ii)+%g-('Q.^.ii)l, ,58) 

_7l V7- 7- 7-/ 7+ V7+ 7+ 7+/. 

where we used a shorthand 7± = 7csc + 7(<I'i ± *I'2) with 7osc and 7(A$) defined in 
Eq. (|49|l . At r = 0, $1,2 — 0, and in the absence of time-dependent perturbations 

= 1. Below we discuss the properties of the functions in various regimes. 

In the limit of high temperature, T ^ 7esc, we obtain 

Q±(a:,y,z)«^-=1=. (59) 

The equality between functions means that the conductance is symmetric with 
respect to magnetic field inversion. However, in low temperature limit T <^ 7±, we 



Sg^Jg^., - 2 



obtain for strong perturbation Ci ^ 7csc 

Q+(a;,0,z)«-^, Q-(a:,0,z)« ^. (60) 



2V2x ^ ' ' ' 2x 

Equation (|58|l with given by Eq. 16U|I shows an important signature of the 
effect of time-dependent perturbations on the conductance — the violation of the 
Onsager symmetry: 



- 7(2$) > 7csc. (61) 



5g^5g^ V C\ 

This breakdown of the Onsager relation is a simple manifestation of lifting of the time 
reversal symmetry in the system with time dependent Hamiltonian. 
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In the limit of low frequency uj <^ 7cscj the conductance g can be represented as 
the result of averaging of the conductance g{{^i}) at stationary perturbation ipi over 
one period 2tt/uj: 

g^j^ 9{W^m^, (62) 

where g{{^i}) can be calculated according to Eq. H15|) with the scattering matrix 
defined by Eqs. l(T^ and ((T^ at fixed values (pi. Because g{ip) has magnetic field 
symmetry |22|. g is also symmetric with respect to inversion of magnetic field. 
Calculations of Eq. (|55|) at cj ^ 7esc give 

n±i ^^ r '^^'^^ T ( , l ^ exp(-(l + 4xsin^ e/2 sin\/2)0 

Q {x,y,0)= -F^ (^/7t) T— ; — . 2./0 ■ 2./0 (63) 

Jq 'iTT^ 7o 1 + 4a;sm^/2sm C/2 

This expression in the limit of high temperature T 3> 7csc has the asymptote 



TT 



Q^{x,y,0) = ^K{-4x), (64) 

and at zero temperature Q^{x, 0, 0) is given by 
n±f nn^ I E {-Ax) + {1 + 4x)K{-ix) 

Q^{x, 0, 0) = — — , 65 

TT 1 + 4a: 

where K{x) and E{x) are the elliptic integrals of the first and second kind respectively. 

Suppression of the conductance fluctuations by slow field u) <C 7esc, see Eq. 
and is the consequence of averaging of the stationary conductance g{(p) over 
different configurations of the full Hamiltonian along the closed contour in parameter 
space, see Fig. 1 b). Thus, the observed d.c. conductance g, Eq. (|^ . is already 
partially averaged over ensemble of random Hamiltonians H and its fiuctuations 
are reduced. As the strength of the perturbation Ci increases, more statistically 
independent configurations of Hamiltonian H contribute to the conductance g and 
fluctuations of g become suppressed. 

However, low-frequency perturbations do not affect the weak localization 
correction to the conductance Ag^i, Eq. H50I) . which is deflned as the difference between 
the averages over orthogonal and unitary ensembles. Only perturbations at frequencies 
w ^ 7csc could suppress Ag^i, see e.g. Eq. (IK^ . when the conductance g, Eq. (fT^ . is 
no longer related to the stationary conductance g^f). In this case the suppression of 
both conductance fluctuations and the weak localization correction to the conductance 
are qualitatively similar and can be interpreted as dephasing. 



5. Photovoltaic current 

Photovoltaic current averaged over ensemble of random Hamiltonian H is zero, 
because there is no speciflc direction for the current to flow. However, for each 
particular conflguration of the quantum dot, a flnite current can flow in either 
direction. To characterize the value of this current, one can find | 61[I47| 

27T j LO 00 00 

var/p, = ^^ I ^^^ Jjljr I deK+{t,t',r,0)B{t-e,t'-9,r), (66) 

''''0 t/2 
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Figure 5. The diagrams for tlie calculation of the variance of conductance 
fluctuations. 



where is defined in Eq. (|5fi|l and 

oo 

B{t,t\T)^^i,f{T) I d^d^'Vit,t-^,T)Vit',t'~^',r) 



. . Tcsc 



[^^,ir,t-0'fj{r,t-0 




'Tcsc 



(67) 



with Lpi{T,t) introduced in Eq. (|49|l . In Fig. |^ we present only the diagram which 
survives at high temperatures and the fuU set of diagrams contributing to Eq. ll^ can 
be found in Refs. |^07], We emphasize that the photovoltaic current has no symmetry 
with respect to inversion of magnetic field, /ph($)/ph(— = 0, this statement in the 
diagrammatic language means that there is no counterpart of the diagram in Fig. (HJ 
that contains the Cooperons; cf. Fig. [3 

Function B{t,t',T), Eq. is related to the electron distribution function in 

the dot. Particularly, at high temperature T ^ Th for harmonic perturbations at 
frequency uj 

B„,,,.i:«,.,(i±^),,(i±^)(~)\ 

corresponds to a square of the distribution function F^^, Eq. (|26|l . Here we introduced 
temperature scale Th according to 



Tv, — wmax 




(69) 



For a single perturbation at high frequency u ^ 7csc with strength Ci and at low 
temperature T <C Th, we can estimate the integrals over ^ and ^' in Eq. H67() as 



/•OC 

/ V{t,t-^,2T)d( 
Jo 



^GSC 



2 — ' (70) 
/o 7csc + 2Cism ujT 

and the function B{t, t' , r) acquires extra factors Eq. (|70|l . which vanish at r ^ 1/Th. 
Narrowing of the distribution functions in time representation means the broadening 
of electron distribution in the energy space. Indeed, energy of an electron in the dot 
changes due to the external field. Such changes result in the redistribution of the 
electrons in the energy space and the new distribution function becomes wider than 
that of electrons in the leads at temperature T. 



Quantum chaotic scattering in time- dependent external fields 



19 



A 



R 



R 



A 



R 



A 



<> <> 



A 



R 



^ = f{t-f)z^v,Mt)-^,{t')] 



Figure 6. The diagram representing the contribution to the variance of the 
photovoltaic current, var/pi^, at high temperature T. 



The new width of the electron distribution function can be estimated from the 
following argument. The transitions occur with rate C\, and electron stays in the dot 
for l/7csc time, so that it experiences of the order C'l/jcsc transitions. After each 
transition electron energy changes by ztcj, and assuming that its motion in the energy 
space can be described by a random walk, we find that on average electron energy 

which is consistent with the estimate Eq. H7(J|I . The 



changes by Th 



above estimate of Th has a meaning only for strong fields, C ^ 7esc) so that the 
diffusion picture in energy space is valid. Otherwise, electrons experience only a few 
transitions with energy change to. Note that function in Eq. H7()(l is periodic in r with 
period ~ 1/w, i-C. at longer times the electron diffusion is no longer described by 
random walk and some structure in the distribution function appears [lOD] . 

For weak harmonic perturbation Cy <ti 7osc, we expand to first order in Cy 
and neglect the second term in Eq. (|67|l . since Th <C T. As the result we obtain 
for the harmonic perturbation, characterized by two functions (pi(t) = cos cut and 
tp2{t) — cos(cjt + 4>): 

oo +e 

C'{{2uj9 - sin 2Lue) + sin 2^6* 



var/ph 



2 2 

e to 



IcscdOe 



-27cBce 



drBir)- 



(71) 



where the linear, Ci, and circular, Cc, perturbation amplitudes were introduced 
according to 



Ci = Cii + 2Ci2 cos(j) + C22, Cc = 2sm(j)JCiiC22 - 



(72) 



Note that the separation of the variance of the current into circular and linear 
contributions corresponds to the classification of the current components introduced 
in Ref. |f)5j . where current through a microjunction in microwave field with linear and 
circular polarizations was studied. In the case of temperature T larger than the escape 
rate, T ^ 7csr-, we have 



/CSC ; 

^2, ,2 



var/nh = 



TT e^uj^ NiN, 7o 



1 



12 47r2 Nl T 



^2 , 

/CSC 



^2 

/CSC 



(73) 



The second term of Ea. (|73|) survives the limit w — > 0, thus reproducing the known 
result for adiabatic pumping |57[ . On the other hand, this term vanishes at high 
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frequency. The term is quadratic in frequency at small frequency and tends to a 
constant at large frequency. 

The linear pumping amplitude Ci in the case of two pumps has the form of 
Ea. H72|) . which implies that the amplitude C\ is just a vector sum of different pumps 
in the parameter space. On the other hand the circular amplitude is determined by 
"uncorrelated" components of matrices Vi,2 and vanishes if Vz (x Vi. 

To describe the variance of the photovoltaic current, we first consider the adiabatic 
limit, when parameters Cij have a special form Cn = C22 — C and C12 = 0. In this 
case the expression for the variance can be written in the compact form for T ^ 7esc 
e^w^ ^jjv^ ^^^^ 2C + (70SC - V 'Tcsc 

As temperature drops down to T ^ 7esc — Nch^i/^n, the variance of /ph saturates to 
""'^'^ = ^^V7e.c(7l + 4C)3- 

The authors of Ref. (^21 showed that at strong perturbation the variance of the 
photovoltaic current in the adiabatic limit is proportional to the perimeter of the 
contour integral in the parameter space Fig. 1 b); var/ph oc VC. This perimeter 
law is the consequence of the lack of correlation between distant points of the 
contour in the parameter space. The total contribution to the pumped current 
consists of uncorrelated contributions of the loop and is proportional to the number 
of independent contributions -y/ C /jcsc ■ In the opposite case of weak perturbation 
C <C 7esc the current /ph is determined by Eq. , see [HII and is proportional to Cc 
(var/ph oc C^). Equations H74|) and (|75|) are consistent with the above arguments for 
power dependence of var/ph. 

When only one perturbation (pi(t) with power Ci is applied, the photovoltaic 
current vanishes in the adiabatic limit. In this case the photovoltaic current is 
quadratic in frequency uj for lu <ti 7csc- For weak pumping var/ph is determined 
by Eq. 1)73(1 with Cc = for arbitrary frequency u. For strong pumping Ci ^ 7osc, 
but still low frequency limit w^Ci <C ^^^c have 

25 e^u^NiN, ^^^^ / c\ 

In the limit of high frequencies, T ^ lo ^ 7csc, the variance of the photovoltaic current 
is given by 

var/r,h — ^ , . (77) 

' 24^ Nl^ T V7csc(7csc + 2Ci) ^ ' 



In the limit of strong pumping this expression has the vCT asymptotic behavior. 

The results for var/ph at finite frequency have the following interpretation. Based 
on Eq. H21|l. we can represent the photovoltaic current in the form similar to Eq. H22|l . 
where the contour of integration is considered in the phase space. The phase space 
includes both parameters ^i{t) and their time derivatives as follows from Eq. 
At weak perturbation with a single parameter (pi{t) = coswt, the contour is an ellipse 
with semi-axes proportional to \/C\ and uj\fC\\ then /ph oc (e.u))ujC\ and var/ph is 
consistent with Eq. 1(73(1 at <C 7osc- In the limit of strong perturbation at low 
frequency w^Ci <C 7esc? the contour in phase plane is long along the 1^9 1 axis but narrow 
in the ^\ direction. The variance of the photovoltaic current is determined by a sum 
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of independent contributions from pairs of the contour along the ip axis, the number 
of these pairs can be estimated as \/C\p^fesc- Each pair consists of two adjacent pieces 
of the contour shifted with respect to each other along (^i [t) axis and contributes as 
ujy/C\ to the total current. As a result, we obtain var/ph oc (ew)^(jj^Cf Finally, 
if the amplitude of the field C or the frequency uo increases further, liPC\ ^ lcsc^ 
the contour does not have adjacent parts and each part of the contour gives an 
independent contribution. Since the number of this parts is C1/70SC1 the variance of 
the photovoltaic current is proportional to (ecj)^\/C', see Eq. (I77II . 

As frequency lo and power increase further, the heating effects become important. 
At Th ^ T, the variance of the photovoltaic current can be roughly estimated if 
electron temperature T in the leads is replaced by T^. Particularly, from Eq. ()77|l we 
obtain the characteristic scale for var/ph: 

-r/ph^^^ — . (78) 

A numerical analysis |83| shows that in fact at C\ ^ 7osc7 the variance of the 
photovoltaic current on the perturbation power has a very weak (log-like) dependence 
on Ci, with typical of value var/ph consistent with the estimate of Eq. H78|) . 

The heating effects manifest themselves in the noise of the photovoltaic current 
as well |77| . In the limit of strong perturbation C\ 3> 7csc at high frequency oj S> 7osc 
the ensemble averaged value of S-p, Eq. (|^ . is 

Sp (X ffciTh- (79) 

The noise of the photovoltaic current has a form similar to the expression for 
the Nyquist- Johnson noise, see Eq. H31|) : the current noise correlation function is 
determined by the conductance of the dot ^ci, and the effective electron temperature. 
Due to the heating by a strong perturbation, the electron distribution function is 
broadened and the new energy scale for the electron distribution function is given 
by Th, see Eq. (f7^ . Thus, the noise of the photovoltaic current averaged over the 
ensemble has a similar origin with the Nyquist- Johnson noise and is determined by 
thermal fluctuations of electrons in the dot out of equilibrium. 



6. Conclusions 



In summary, we reviewed the random matrix description of electron transport through 
an open quantum dot, subject to time-dependent perturbations. We expressed the dc 
current through the dot in terms of the scattering matrices, and considered such 
components of the current as the photovoltaic current, independent from the bias 
voltage, and the linear in the bias current, characterized by the conductance. The 
scattering matrices are calculated in terms of time-dependent Hamiltonian, that 
belongs to a Gaussian ensemble of random matrices. We then presented the diagram 
technique to perform ensemble averaging and applied this technique to calculate 
different statistical properties of the electron transport through the dot. 

The main results can be summarized as follows. The weak localization 
correction to the conductance and conductance fluctuations are both suppressed 
by time-dependent perturbation. However, the suppression has different parametric 
dependence on perturbation frequency. The photovoltaic current can be represented as 
a sum of circular and linear terms. These term have different frequency dependence: 
the circular term dominates at low frequencies and represents the adiabatic charge 



Quantum chaotic scattering in time- dependent external fields 



22 



pumping, while the linear term dominates at high frequencies. The photovoltaic 
current and its noise are determined by the actual width of the electron distribution 
function in the dot, on the other hand, the variance of the conductance fluctuations 
is determined by electron temperature in the leads. These results are in qualitative 
agreement with experiments, described in Refs. |10[ I82| . 

We described calculations using the Hamiltonian approach to the statistical 
description of the electron transport, a detailed description of the scattering matrix 
approach for time-dependent system can be found in Ref . 47„ , where the same results 
were obtained. 

In this review we considered the electron system neglecting the interaction 
effects and assumed spin degeneracy. The effect of electron-electron interaction 
can be disregarded only in the limit of the large number of open channels, but 
as the number of open channels decreases, the interaction effects become more 
important |84[ 1861 1871 I1(J1| . The interplay of the interaction and time-dependent 
perturbation was addressed in Refs. |1U2[ ll(J3| . 

In semiconductor quantum dots in absence of magnetic field electron spin states 
are nearly degenerate. However, if magnetic field is applied, the spin degeneracy is 
lifted and currents of electrons with opposite spin orientations are not identical. In this 
case a spin current can be generated by time-dependent perturbation, similar to the 
photovoltaic charge current |ll)4| . this effect was studied experimentally in Ref. |ll)5| . 
Another modification of the system, considered in the present review, is a quantum 
dot connected to superconducting leads and was studied theoretically in 106.. .lOTj , 
the experimental realization of such a system remains a challenging task. 
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Appendix A. 

We denote the wave function of electrons in channel a by tpaixjt) with a; < for 
incoming electrons and x > for outgoing electrons, see Fig. 1 a). The boundary 
a; = is described by a superposition of the incoming and outgoing electron states 
and we denote it by ipa{0, t)- The wave function of electrons in the dot is denoted by 



which is defined in terms of the retarded, advanced and Keldysh components as 



We introduce the matrix Green function 




(A.l) 



gl^f^{t,t',x,x') 

gifit,t',x,x') 



~zeit-t'){{^P„{x,ty,i;l{x',t')}) 
ieit'-t){{i;^ix,t);i^lix',t')}), 
-t{[i:a{x,t);^Ux\t')]), 
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where [A; B] — AB — BA and {A; B} — AB + BA. The similar expressions can be 
written down for Qia{t,t\x') Green function, with ^pai^it) replaced by ipiit)- 

For non- interacting electrons, moving towards the dot {x; x' < 0), the Green's 
function is: 



Gap 



(f i' T t') — f '^afii^ ^ ^' ~ C!^p{t - t',X - x') \ 

^ ' ' ^ Gi^it-t',x-x') J 



(A.2) 



where 

Gf^(t, x) = te{t) 5^p5 [vpt - x) , (A.3) 

Gipit, x)= - iQ{-t) So^pS {vpt - x) , (A.4) 

G^p{e,x) = Ue) {G^^pie.x) - Gi^ie^x)) , (A.5) 

and /(e) is the distribution function of electrons in channel a. If incoming electrons 
are in equilibrium at temperature T, 

e - eVa 



fa{e) = tanh- 



2T 



(A.6) 



with Va being the voltage applied to the reservoir connected to the dot by channel a. 

+1 



are 
i 



The equations of motion for the Green functions Qap{t,t' , x, x') and Gja{t,t' , x' 

g^pit, t', X, x') = six)Wc.^g,f3it, t', x') + Sit - t')5{x - (a.t) 

Gjait^t'.x') = wlgpo.{t,t',Q,x'). (A.8) 



d__ d_ 

dt di 
d 

^gl-H.,it) 



Due to causality principle, G'^^(t, i', 0, x') = for x' < 0. This observation 
significantly simplifies further calculations. Indeed, we can represent the Keldysh 
component of the Green's function in the left hand side of Ea. (|A.8|) in the form 



(t,t',x')= I dt^ 





1 


1 dti 


id/dt-H{t) 



{t,ti) w]pgo,p{ti,t',o,x') 



(A.9) 



The corresponding advance component is zero. Here \/{id/dt— H(t)) is the retarded 
component of the electron Green's function in the dot. This definition is different 
from that given in the main part of the paper, see Ea. (|ll|l . The latter will appear 
naturally in the end of this calculation with an additional term ~ W^W, see Eq. (|13f) . 
describing escape of electrons from the dot through the leads. We represent Ea. ljA.7p 
in the form 



gly{t,t',x,x') =G%{t~t',x~x')+ dtidt2G^^{t-ti,x) 



W- 



id/dt-H{t) 



(ii,t2)el^^(i2,t',0,x'), 



(A.IO) 



- 7^ 



take the limit = and, using 

G%{t -t',0) from Ea. (fO)l . obtain for a;' < 



gif{t,t',0,x')^ Idt, 



1 - W- 



id/dt~H{t) 



{tM)G%\h.t',Q,x'). 



Substituting this expression to Eq. ljA.10|) and taking x = +\6\ 0, wc find 



QZ\t,t\+\5lx') ^ ldtiSo,-,{tM)G'^p{ti-t',-x'), x'<0, 



(A.ll) 



(A.12) 
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where the scattering matrix S^/sit, t') is defined by Eq. (|12|) . 

Equation ljA.12p is vahd for x' < 0. We have to repeat the procedure described 

above to calculate the electron Green's function in the leads for x' > 0. Since the 

equations which determine evolution of the Green's function from x' < to a;' > are 
conjugated to those for x, we obtain 

Gifit,t',+\S\,+\S\) = JJdtidt2So.^it,ti)G'^siti~t2,0)Slit2,t'). (A.13) 
The currents in the left (right) leads are given by 

(/i(,)(t))=e^;p J2 {sfJit.t,+\6\,+\6\)-gi''Jit^t,-\6\,-\6\)), (A.14) 

where a — 1 . . . Ni for left lead and a = A^i + 1 . . . iVch for right lead; coordinate 6 
is in the lead just before the contact with the dot: 6 > and (5-^0. Function 
Gaa\t,t, —\S\, —\6\) is taken for incoming electrons and is given by Ea. ljA.5|l and 
consequently, 

/ + 00 7 
e^-'fic)^. (A.15) 

The total dc current from the dot should be zero + {Ir) = to ensure no 
charge accumulation on the dot, where (/i(i-)) = {I\{j^{t))dtlTo. Therefore, we can 
rewrite the expression for dc current through the dot as 

m - . (« . _(,,). ,A,e) 

Substituting Eqs. I|A.13|) and (|A.15|) into Eqs. (jA.14|) and using Eq. (|A.16p we obtain 
Eq. @. 



Appendix B. 

In this Appendix we derive Eq. H29|l for the current noise correlation function through 
a quantum dot. The quantum mechanical operator of the current through left (right) 
lead is 

/i(,)(t) = ei;F {€{t,+6)Mt,+S)-^i{t,-S)Mt,~5)) (B.l) 

with -tpaltjiS) being the operator for outgoing (+6) or incoming (—S) electrons 
through channel a, cf. to Eq. (|A.14p . 

Substituting the expression for the current operator, Eq. (jB.l|) into Eq. H28|l and 
using the charge conservation in the dot on time Tq 3> I/ijJ, wc obtain the following 
expression for the current correlation fimction (below S ^ 0, but i5 > 0): 

S = e^vl ° dtdt' (tr I A^< {t\ t, +5, +5)kg> {t, t', +6, +S) } 

- tr {a^< {t\ t, +5, -S)Ag>{t, t', -S, +,5)} 
^ tr {Ag<(t', t, -S, +S)Ag>{t, t', +S, -5)} 

+ tr {a^< {t\ t, -5, -S)Ag>{t, t\ -5, -5)} ) . (B.2) 
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Here we introduced the electron Green's functions in the leads according to the 
following definitions (for a review of the Keldysh Green's function formalism see f2^)■. 



g<^it, t\ X, x') = iii^lit', x')^„{t, x)), g>p{t, t', X, x') = -iiiPo.it, x')). (B.3) 



The Green functions g^'^ can be written in terms of the retarded, advanced and 
Keldysh Green's functions: 



The next step is to represent the Green's functions as a product of incoming 
electron Green's functions, Eqs. ljA.2p - (|A.5p . and the scattering matrix, Eq. (|12|) . The 



Now the derivation of Eq. (|29ll reduces to simple algebraic calculations. With the 
help of Eqs. (|B.4|I and (|B.5I) we rewrite Eq. (|B.2p in terms of the retarded, advanced 
and Keldysh components of the Green function. Then we represent these components 
as a product of scattering matrices and the Green's functions of the incoming electrons, 
using Eqs. (ESJ-lEIinj- The result is given by Eq. 
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